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Abstract 

The diffractive photoproduction of p mesons, ep — ► epY, with large momentum trans- 
fer squared at the proton vertex, \t\, is studied with the HI detector at HERA using an 
integrated luminosity of 20.1 pb _1 . The photon-proton centre of mass energy spans the 
range 75 < W < 95 GeV, the photon virtuality is restricted to Q 2 < 0.01 GeV 2 and the 
mass My of the proton remnant is below 5 GeV. The t dependence of the cross section 
is measured for the range 1.5 < \t\ < 10.0 GeV 2 and is well described by a power law, 
dcr/d|i| oc \t\~ n . The spin density matrix elements, which provide information on the he- 
licity structure of the interaction, are extracted using measurements of angular distributions 
of the p decay products. The data indicate a violation of s-channel helicity conservation, 
with contributions from both single and double helicity-fiip being observed. The results are 
compared to the predictions of perturbative QCD models. 
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1 Introduction 



Diffractive vector meson production in ep interactions with large negative four-momentum 
transfer squared at the proton vertex, t, provides a powerful means to probe the nature of the 
diffractive exchange. It has been proposed as a process in which Quantum Chromodynamics 
(QCD) effects predicted by the BFKL evolution equation [1,2] could be observed. Theoreti- 
cally, the large momentum transfer provides the hard scale necessary for the application of per- 
turbative QCD (pQCD) models. In such models, diffractive vector meson production is viewed, 
in the proton rest frame, as a sequence of three processes well separated in time: the intermedi- 
ate photon fluctuates into a qq pair; the qq pair is involved in a hard interaction with the proton 
via the exchange of a colour singlet state, and the qq pair recombines to form a vector meson. At 
leading order (LO), the interaction between the qq pair and the proton is represented by the ex- 
change of two gluons. Beyond LO, the exchange of a gluon ladder has to be considered which, 
in the leading logarithm (LL) approximation, can be described by the BFKL equation. The 
treatment of the formation of the vector meson, which is a non-perturbative process, involves a 
parameterisation of the meson wavefunction. 

At high \ t\, the t dependence of vector meson production and of the spin density matrix elements 
has been measured in photoproduction by ZEUS [3] (p, 0, J/ip) and HI [4,5] (J/ip, ip(2s)). HI 
has also measured the t dependence of the spin density matrix elements for the p meson at high 
\t \ in electroproduction [6]. The data on light vector meson production (p, 0) at large \t\ indicate 
a violation of s-channel helicity conservation (SCHC), i.e. the non-conservation of the helicity 
between the exchanged photon and the vector meson. This is in contrast to measurements of 
the heavier J /ip meson, where s-channel helicity is seen to be conserved [3,4]. 

This paper presents new measurements of the diffractive production of p mesons at large \t\, in 
the range 1.5 < |i| < 10 GeV 2 : 

ep — * epY; p — > 7r + 7r~, (1) 

in the photoproduction regime, i.e. Q 2 ~ GeV 2 where Q 2 is the modulus of the squared four 
momentum carried by the intermediate photon. The system Y represents either an elastically 
scattered proton or a low mass dissociated system. A clean signature and low background rates 
make this process experimentally attractive and it is possible to precisely determine the kine- 
matics of the process through an accurate measurement of the vector meson four momentum. 
The dependence on t of the cross section is measured and the spin density matrix elements are 
extracted, which provide information on the helicity structure of the interaction. 



2 Perturbative QCD Models 

Perturbative QCD calculations at leading order (two-gluon exchange) have been performed 
for high \t\ vector meson photoproduction in [7]. For light vector mesons, these calculations 
indicate that, although the initial photon is transversely polarised, the vector meson is produced 
predominantly in a longitudinal state. Intuitively, this can be understood as follows. Since the 
light quark and antiquark have opposite helicities (chiral-even configuration), the pair must be in 
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an orbital momentum state with projection L z — ±1 onto the photon axis, in order to conserve 
the photon spin projection. The hard interaction between the dipole and the gluon system does 
not affect the dipole size nor the quark and antiquark helicities, but modifies the dipole line of 
flight. This implies a damping by a factor oc l/\t\ of the probability of measuring the value 
L z = ±1 for the projection of the dipole angular momentum onto its line of flight, and hence 
of transversely polarised vector meson production. 

However, in contrast to these pQCD expectations, experimental observations indicate that light 
vector mesons are produced predominantly in a transverse polarisation state [3] . As first realised 
in [8] an enhanced production of transversely polarised p mesons can be accounted for by the 
possibility for a real photon to couple to a qq pair in a chiral-odd spin configuration even in the 
case of light quarks 1 . 

The data presented in this paper are compared to two theoretical predictions: a fixed order cal- 
culation in which the hard interaction is approximated by the exchange of two gluons 2 , and a 
LL calculation in which it is described according to the BFKL evolution. In both cases, the 
chiral-odd component of the photon wavefunction is obtained by giving the quarks a "con- 
stituent quark mass" m = my/2, where m v is the mass of the vector meson, and a set of QCD 
light-cone wavefunctions for the vector meson [1 1] is used. 

The BFKL calculation is described in [9, 10]. The BFKL resummation cures some possible 
instabilities which might affect the two-gluon predictions [1]. The expansions on the light-cone 
of the relevant hadronic matrix elements are performed up to twist-3, i.e. next-to-leading twist. 
The leading logarithm nature of the calculation prevents an absolute prediction for the normal- 
isation of the cross sections, due to the presence of an undefined energy scale A. In [9, 10] this 
scale is allowed to run with t according to A 2 = rriy — ^t, where 7 is a free parameter. The value 
of the strong coupling constant is fixed since this is appropriate in the LL approximation and 
has proved successful [12, 13] in describing previous data. The parameter values, as obtained 
from a fit to the ZEUS data in [3], are: 7 = 1.0 and af FKL = 0.20 [14]. 



3 Data Analysis 
3.1 Event Selection 

The data used for the present analysis were taken with the HI detector in the year 2000 and 
correspond to an integrated luminosity of 20.1 pb _1 . In this period, the energies of the HERA 

'Note that in the case of heavy vector mesons, the relevant non-relativistic wavefunction, with equal sharing of 
the photon longitudinal momentum by the two heavy quarks, ensures that the amplitude for producing a longitu- 
dinally polarised meson vanishes. Only the chiral-odd amplitude, which is naturally present due to the non-zero 
quark mass, remains. Hence heavy vector mesons are expected to be transversely polarised, as confirmed by the 
experimental observations. 

2 The two-gluon model predictions for the kinematic region considered here were provided by the authors 
of [9, 10]. The model differs from that proposed in [8] in the way the chiral-odd configurations are introduced. 
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proton and positron beams were 920 GeV and 27.5 GeV, respectively. The kinematic domain 
of the measurement is: 



Q 2 < 0.01 GeV 2 

75 < W < 95 GeV 

1.5 < \t\ < 10 GeV 2 

M Y <5GeV (2) 

where W is the photon-proton centre of mass energy and My is the mass of the proton remnant 
system. 

The relevant parts of the detector, for which more details can be found in [15], are the central 
tracking detector (which covers the polar angular range 20° < 9 < 160°), the liquid argon 
(LAr) calorimeter (which covers the polar angular range 4° < 9 < 154°) and an electron tagger 
located at 44 m from the interaction point, which detects the scattered positron at a small angle 
to the backward direction 3 . The 44 m electron tagger is a 2 x 3 array of Cerenkov crystal 
calorimeters used to select photoproduction events in the range Q 2 < 0.01 GeV 2 . The trigger 
system used in this analysis selects events with an energy deposit greater than 10 GeV in the 
44 m electron tagger, at least one charged track with a transverse momentum above 400 MeV 
in the central tracker and a reconstructed interaction vertex. Additionally, there is a veto on the 
amount of energy deposited in the forward region of the LAr. 

Events corresponding to reaction (1), in the kinematic range defined by relations (2), are se- 
lected by requesting: 



• the reconstruction of an energy deposit of more than 15 GeV in the 44 m electron tagger 
(the scattered positron candidate); 

• the reconstruction in the central tracking detector of the trajectories of exactly two oppo- 
sitely charged particles (pion candidates) with transverse momenta larger than 150 MeV 
and polar angles within 20° < 9 < 155°. In order to ensure a well understood trigger 
efficiency, at least one track is required to have a transverse momentum above 450 MeV; 

• the absence of any localised energy deposit larger than 400 MeV in the LAr calorimeter 
which is not associated with either of the two reconstructed tracks. This cut reduces 
backgrounds due to the diffractive production of systems decaying into two charged and 
additional neutral particles. Further, it limits the mass of the proton dissociative system to 
My < 5 GeV which ensures that the events lie within the diffractive regime My <C W. 
Restricting the analysis to low values of My also reduces the uncertainties arising from 
the parametrisation of the My dependence of the cross section; 

• events in the mass range 0.6 < M nn < 1.1 GeV and discarding events with 
Mkk < 1-04 GeV, where M nw and Mkk are the invariant masses of the two selected 
tracks when considered as pions or kaons respectively (no explicit hadron identification 
is performed for this analysis). The latter cut reduces the background due to diffractive 
production of <fi mesons. 

3 In the HI convention, the z axis is defined by the colliding beams, the forward direction being that of the 
outgoing proton beam and the backward direction that of the positron beam. Transverse and longitudinal momenta 
are defined with respect to the proton beam direction. 
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The final sample consists of 2628 events. Further details of this analysis may be found in [16]. 



3.2 Kinematics and Helicity Structure 

The three momentum of the p meson is computed as the sum of the two charged pion candidate 
momenta. The variable W is reconstructed from the p meson, rather than from the energy of 
the scattered positron candidate which is less precisely measured by the electron tagger, using 
the Jacquet-Blondel method [17]: 

W 2 ~2E p (E p -p z , p ) (3) 

where E p is the energy of the incoming proton and E p and p Z:P are the energy and longitudinal 
momentum of the p meson, respectively. In the photoproduction regime, the variable t is well 
approximated by the negative transverse momentum squared of the p meson, t ~ — p\ p . 

The measurement of the production and decay angular distributions provides information on 
the helicity structure of the interaction. Three angles are defined [18] as illustrated in Fig. 1: $ 
is the angle between the p production plane (defined as the plane containing the virtual photon 
and the p meson) and the positron scattering plane in the ^p centre of mass system; 9* is the 
polar angle of the positively charged decay pion in the p rest frame with respect to the meson 
direction as defined in the ^p centre of mass frame and 0* is its azimuthal angle relative to the 
p production plane. 

In this paper, the distributions of the angles 0* and 9* are analysed (the angle $ is not acces- 
sible in photoproduction) giving access to the spin density matrix elements r^, Re [r^ 4 ] and 
[19]. For a p meson decaying into two pions, the normalised two-dimensional angular 
distribution [19], averaged over $, is 



d 2 a 3 



^(l-0 + ^(3r o 4 -l)cos 2 r 



erdcos0*d0* An 

-V2Re [rft] sin 29* cos 0* - r? 4 : sin 2 9* cos 20* 



(4) 



Integrating over cos 9* or 0* further reduces this distribution to the one dimensional distributions 

da 



and 



d cos 9* 

da 



ocl-r ° 4 + (3r ° 4 -l)cos 2 r (5) 



oc 1 -2r° 4 -, cos 20*. (6) 
d0* 1-1 

The spin density matrix elements are defined as bilinear combinations of the helicity amplitudes 
M\ \ v , where A 7 , \y = — , 0, + are the respective helicities of the photon and the vector me- 
son. For photoproduction, where the photon is quasi-real, the longitudinal photon polarisation 
component is negligible and only the transverse polarisation states remain. The photon-meson 
transitions can thus be described in terms of three independent helicity amplitudes M ++ , M +0 , 
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M + _ 4 , which correspond to no change in helicity (no-flip), a single change in helicity (single- 
flip) and a double change in helicity (double-flip), respectively. In this case, the matrix elements 
are related to the helicity amplitudes by 

r 04 = l M +o! 2 (7) 

00 " |M ++ | 2 + |M +0 | 2 + |M+_| 2 V ; 

i m ++ m; - m + _m; 

10 2|M ++ | 2 + |M +0 | 2 + |M + _| 2 V ; 



1 M ++ M* + _ + M + _M* ++ 
1-1 2|M ++ | 2 + |M +0 | 2 + |M 4 I 2 ' 



Under s-channel helicity conservation, only the amplitude M ++ is non-zero and consequently 
the Tqq, ?~io and r°i 1 matrix elements should all be zero. In contrast, both the BFKL and two- 
gluon models predict a violation of SCHC. In the case of the LL BFKL model, the helicity 
amplitudes are predicted to follow a hierarchical structure with |M ++ | > |M+_| > |M +0 | [9, 
10]. 



3.3 Monte Carlo Simulation 

A Monte Carlo simulation based on the DIFFVM program [20] is used to describe the diffrac- 
tive production and decay of p mesons, and to correct the data for acceptance, efficiency and 
smearing effects. The Q 2 and W dependences of the cross section are taken from previous mea- 
surements [21]. The t dependence is taken according to the power law measured in the present 
analysis following an iterative procedure. The simulation includes the angular distributions 
corresponding to the measurements of the present analysis for the r^, r^_i and the Re [r^] 
matrix elements. The M Y spectrum is parameterised as da/dMy oc /(My)/My 15 [22], where 
/(My) = 1 for My > 3.6 GeV 2 and, at lower masses, is a function which accounts for the 
production of excited nucleon states. Other angular distributions and correlations are taken in 
the s-channel helicity conservation approximation. The mass distribution is described by a rel- 
ativistic Breit-Wigner distribution, the mass and the width of the p being fixed [23], including 
skewing effects resulting from the interference with open pion pair production [24, 25] taken 
from the current analysis. For studies of systematics uncertainties, all simulation parameters 
have been varied within errors (see section 3.5). 

The DIFFVM simulation is also used for the description of the u, 4> and p' backgrounds (see 
next section). Here the t distributions are described using the same parameterisation as for the 
p meson and the angular distributions are kept in the SCHC approximation, which is justified 
by the smallness of these contributions. 

All generated samples are passed through a detailed simulation of the detector response based 
on the GEANT3 program [26], and through the same reconstruction software as used for the 
data. Figure 2 presents the observed and simulated distributions for several variables of the se- 
lected sample of events. The simulated distributions, which include the amounts of background 

4 The three corresponding amplitudes M , M_o and M |_ are not independent since they satisfy M ++ = 

M , M + _ = M-+ and M +0 = -M_ due to parity symmetry. 
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discussed in section 3.4, are normalised to the number of observed events. Reasonable agree- 
ment is observed for all distributions, showing that the Monte-Carlo simulations can reliably be 
used to correct the data for acceptance and smearing effects. The structure in the 0* distribution 
(Fig. 2f) is due to the low geometrical acceptance of the central tracking detector in the case 
where the angle between the p meson production and decay planes is small ((f)* ~ 0° or 180°), 
which leads to the emission of one of the pions at a small angle relative to the beam direction. 



3.4 Backgrounds 

Diffractive photoproduction of uo, and p' mesons 5 can fake p production through the decay 
channels: 

uo — > 7T + 7r~7r°, 

p'^p^V, p ± ^7T ± 7T , (10) 

if the decay photons of the n° or the K\ mesons are not detected 6 . This happens in the cases 
where the energy of the neutral particle is deposited in an inactive region of the detector, is 
associated to the charged pion tracks, or is below the noise threshold. Diffractive photoproduc- 
tion of uo and mesons also gives the same topology as that of the p meson within the detector 
through the decay channels 

uo — > 7T + 7r~. 

-> K + K~. (11) 



To estimate the corresponding backgrounds, the uo, and p' cross sections were taken from 
measured ratios to the p cross section in the Q 2 range relevant for the analysis: a u / a p — 
0.106 ± 0.019 [27], a 4, / a p = 0.156±g;g?| [3] and a(p' -> pVir ) / a p = 0.2 ± 0.1. In the 
latter case, the ratio is obtained from the ratio a(p' — > p°n + n~) / a p — 0.10 ± 0.05, measured 
in electron [28] and muon [29] scattering off a liquid hydrogen target, under the assumption 

a(p>^p+n^) + a(p>^p-n+n°) _^ ^ 
a(p' — > p°7r+7r~) 



The background contributions in the selected kinematic domain (2) and for the selected p mass 
range are estimated and subtracted separately for each measurement interval using the Monte 
Carlo simulations. In total, the backgrounds amount to 0.5%, 0.2% and 1.2% for uo, and p' 
production, respectively. 

5 The detailed structure [23] of the p' state is not relevant for the present study. The name p' is, therefore, used 
to represent both the p'(1450) and the //(1700). In the DIFFVM simulation, the p' mass and width are taken as 
1450 MeV and 300 MeV, respectively. 

6 The contribution from background processes leading to more than two charged particles in the final state is 
negligible. 
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3.5 Systematic Uncertainties 



The systematic uncertainties on the measurements are estimated by varying the event selection, 
the parameters of the p Monte Carlo simulation and the properties of subtracted backgrounds. 
The following sources of systematic error are taken into account: 

• Uncertainties in the p simulation 

The uncertainty on the input t distribution, used to compute acceptances and smearing ef- 
fects and to adjust the measurements to the mean t value for each t interval, is taken into 
account by varying the exponent of the Monte Carlo power law by ±0.5. The uncertainty 
in the modelling of the dissociative proton system Y is estimated by reweighting the pro- 
ton remnant mass distribution by factors (1/My) ±0 3 [22]. For the angular distributions, 
the spin density matrix elements are varied around the values measured in the current 
data according to the spread of the observed results with t by ± 0.03 for r^, ± 0.02 for 
Re [rJo] and for r^. Finally, in the low \t\ region where non-zero skewing is ob- 
served in the p line shape, the skewing parameter is varied according to the uncertainty 
of the fit to the invariant mass distribution (see section 4). 

• Uncertainties on the background distributions 

The amount of background is varied by changing the ratio of the background cross sec- 
tions to the p cross section according to their uncertainties, as quoted in section 3.4. The 
t dependence of the p' distribution, which provides the largest background contribution, 
is further varied using weighting factors of (l/|t|) ±2 . 

• Uncertainties in the detector description 

Uncertainties on the detailed p t and angular dependences of the trigger efficiencies are 
taken into account by varying them within their estimated errors. The uncertainty on the 
tracking acceptance at large angles is estimated by varying the cut on the polar angle of 
the reconstructed tracks between 150° and 160°. The threshold for the detection of energy 
deposits in the LAr that are not associated to the two pion candidates is varied between 
300 MeV and 500 MeV. Finally, the influence of the uncertainty in the description of the 
electron tagger acceptance is estimated by shifting the acceptance range in W by 3%. 

For the measurement of the t dependence, the largest sources of systematic uncertainty are the 
slope of the t distribution in the MC, the variation of the LAr energy threshold and the variation 
of the upper 9 cut. Additionally, for the measurement of the spin density matrix elements, 
the parameterisation of the matrix elements in the MC provides a significant effect. For each 
contribution, the relative effect on the extracted t slope is less than ±1%, and the effect on the 
measured spin density matrix elements is less than ±0.015. 

The total systematic error on the cross section is obtained by adding the individual contribu- 
tions, which are considered as uncorrelated, in quadrature. Correlated systematics which affect 
only the normalisation cancel since only normalised cross sections are presented here. For the 
extraction of the t slope and the spin density matrix elements, the systematic error is obtained 
by repeating the appropriate fit after shifting the data points according to each individual sys- 
tematic uncertainty. Again, the individual contributions are added in quadrature. 
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4 Results 



In each bin of the kinematic variables, the cross section is computed from the number of events 
in the bin, fully corrected for backgrounds, acceptance and smearing effects using the Monte 
Carlo simulations described above. At low \t\, the data are further corrected for the skewing 
effect obtained from a fit to the present data using the Ross-Stodolsky [25] parameterisation 
in which the mass and the width of the p are fixed [23]. This correction amounts to 2.6% 
for 1.5 < \t\ < 2.2 GeV 2 . At higher \t\, the skewing effect is negligible. All cross sections 
presented below, as well as the theoretical model predictions, are normalised to their integrals 
in the respective kinematic domain. 

4.1 Dependence on t 

The t dependence of the ep — > epY cross section is presented in Fig. 3 and table 1. The data 
are plotted at the mean value in each t interval determined according to the parameterisation 
of the t dependence. The data are well described over the measured range of t by a power law 
dependence of the form dcr/d|t| oc \t\~ n where n = 4.26 ± 0.06 (stat.) ( s Y s t-)> as 
determined by a x 2 minimisation. An exponential parameterisation of the form da/d|t| oc e~ 6 '*l 
is unable to describe the data over the full t range. The data in Fig. 3 are compared with the 
predictions of the two-gluon model both with fixed and running a s and with those of the BFKL 
model [9, 10]. The BFKL model provides a reasonable description of the t dependence, in 
contrast to the two-gluon model predictions. 

The ZEUS Collaboration has previously published data on the diffractive photoproduction of p 
mesons with proton dissociation in the range 1.1 < \t\ < 10.0 GeV 2 [3] and observes a power 
law with exponent n = 3.21 ± 0.04 (stat.) ±0.15 (syst.). This is significantly shallower than 
the result obtained here. The difference in slope can be understood in terms of the difference 
in the kinematic region over which the two measurements are made, in particular the maxi- 
mum value of M Y : the phase space of the ZEUS measurement corresponds to M Y < 12 GeV 
(31 GeV) at \t\ = 1.5 GeV 2 (\t\ = 10 GeV 2 ) whereas the phase space of this measurement 
is My < 5 GeV. Both measurements are well described by the BFKL model using similar 
parameters [10, 14]. 

4.2 Spin Density Matrix Elements 

The spin density matrix elements are extracted by a two-dimensional likelihood fit of equa- 
tion (4) to the data. The normalised single differential distributions in cos 9* and <p* are shown 
in Fig. 4 for three ranges of t. The solid curves show the projection of the two-dimensional fit 
and the dashed curves show the expectation of s-channel helicity conservation. A fiat <fi* be- 
haviour is clearly disfavoured, indicating a violation of SCHC. The values of the three extracted 
matrix elements are shown in Fig. 5 and table 2 as a function of \t\. Within the experimental 
uncertainty, no strong dependence on t is observed. Measurements of the spin density matrix 
elements for the photoproduction of p mesons obtained by the ZEUS Collaboration [3] are also 
shown. There is a reasonable agreement between the results of the two experiments. 
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\t\ range 
(GeV 2 ) 


(1*1) 
(GeV 2 ) 


I/a d<7/d|f| 
(GeV" 2 ) 




1 79 


1 I 7« _|_ n noo +0.020 
J — l / u zn u.uoi — 012 




9 9^ 


n AA7 + (1(11(1 +0.0051 

U.^^f in U.UlcJ —0 018 


z.o — o.U 


Z. < O 


n 1 sq 4- n m +0.0072 

U.lOd ± U.U1J _o.0062 


o.U — 0.0 


o./o 


n ns^n -i- n nn7R +0.0025 

U.UoOU it U.UU 1 —o 0041 


3.5-4.0 


3.74 


0.0401 ± 0.0051 +Hmq 


4.0-5.0 


4.45 


0.0188 ± 0.0025 ±8:0000 


5.0-6.0 


5.46 


0.00848 ± 0.00167 t° S 


6.0 - 10.0 


7.56 


0.00185 ±0.00039 IHoooo 



Table 1: The normalised differential cross section for ep — > epY as a function of t. The first 
errors are statistical and the second are systematic. The kinematic range of the measurement is 
given in (2). 

The small values of r^, which is directly proportional to the square of the single-flip helicity 
amplitude M +0 , signify that the probability of producing a longitudinally polarised p from a 
transversely polarised photon is low, varying from (4 ± 2)% at \t\ = 1.79 GeV 2 to (6 ± 6)% 
at \t\ = 4.69 GeV 2 . The non-zero values of Re [r^] confirm that, although small, a single- 
flip contribution is present. The production of transversely polarised p mesons must, therefore, 
dominate and the finite negative values of r ^ show clear evidence for a helicity double-flip 
contribution. Both these observations indicate a violation of the SCHC hypothesis. This is in 
contrast with the results on the J/ip meson, where, within experimental errors, the measured 
spin density matrix elements [3,4] are all compatible with zero. 

The two-gluon model predictions (dotted lines in Fig. 5) [9, 10] are unable to describe the 
measured spin density matrix elements. In particular, the model predicts too high values of 
Tog, i.e. too high probabilities for producing longitudinally polarised p mesons. For the BFKL 
predictions [9, 10], which are shown by the full lines in Fig. 5, is well described but the 
prediction for r 04 i 1 is too negative and the wrong sign for Re [r^] is predicted. The inability to 
describe the Re [r^] matrix element is the major obstacle for the BFKL model. 

5 Summary 

The diffractive photoproduction of p mesons, ep — > epY, has been studied using the HI detector 
at HERA in the kinematic range Q 2 < 0.01 GeV 2 , 75 < W < 95 GeV, 1.5 < |i| < 10 GeV 2 
and My < 5 GeV. The t dependence of the cross section is measured and fitted with a power 
lawoftheform \t\~ n , which fits the data well and results inn = 4.26±0.06 (stat.) (syst.). 
It is reasonably described by a BFKL-based model, while both two-gluon predictions consid- 
ered here, with different treatments of the strong coupling, fail to describe the data. 
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\t\ range 
(GeV 2 ) 


(l*l> 
(GeV 2 ) 


r 04 
'00 


Re [r° 4 ] 


r 04 
'1-1 


1.5 - 2.2 
2.2-3.5 
3.5 - 10.0 


1.79 
2.64 
4.69 


0.038 ± 0.017 1°;^ 
0.029 ± 0.025 t °on 
0.062 ± 0.058 ±g:8ii 


o.o64 ± 0.012 ±8:S?i 
o.o3i ± o.oi9 ±8:g?I 

0.057 ± 0.034 1°;°$ 


-0.088 ± 0.015 toml 
-0.138 ±0.021 t °on 

-0.H9 ±0.044 ±g : g55 



Table 2: The three spin density matrix elements 7q 4 , and Re [r^ 4 ] for p meson photopro- 
duction as a function of \t\. The first errors are statistical and the second are systematic. 



The spin density matrix elements Tqq, r^_ 1 and Re [rj 4 ] are measured as a function of t. The 
r^ti and Re [r^ 4 ] matrix elements differ significantly from zero, thus confirming the violation 
of s-channel helicity conservation, with contributions from both single and double helicity-flip 
observed. The models considered are unable to describe the spin density matrix elements. The 
two-gluon model predicts far too large a probability of producing a longitudinally polarised p 
meson, given by the matrix element. While the BFKL based model is able to describe the 
matrix element well, the prediction for r° 4 : is too negative and Re [r^ 4 ] has the wrong sign. 
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Rest Frame 

Figure 1 : Production and decay angles used to analyse the polarisation of the p meson. 
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Figure 2: Distributions of the selected events in the kinematic domain (2) and the invariant 
mass range 0.6 < < 1.1 GeV: a) transverse momentum of the pions; b) invariant mass of 
the two pions; c) modulus of the square of the four momentum transferred at the proton vertex 
t; d) photon-proton centre of mass energy W; e) cosine of the polar angle 9* of the positively 
charged decay pion in the p rest frame and f) its azimuthal angle <p* the •yp centre of mass frame. 
The points represent the data and the histograms show the Monte-Carlo predictions normalised 
to the data, including the u, 4> and p' backgrounds (filled histograms). 
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|t| (GeV 2 ) 

Figure 3: The t dependence of the ep — » epY cross section. The inner error bars show the 
dominating statistical errors, while the outer error bars represent the sum of the statistical and 
systematic errors added in quadrature. The dashed line is the result of a fit to a power law 
distribution \t\~ n , which results in a power n = 4.26 ± 0.06 (stat.) 1^4 (syst.). The full 
line shows the prediction from the BFKL model and the dotted and dashed-dotted lines show 
the predictions from the two-gluon model with fixed and running a s , respectively. 
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Figure 4: Normalised decay angular distributions for p meson photoproduction in three bins 
of \t\. The left column (a,c,e) shows the polar distribution cos#* and the right column (b,d,f) 
shows the azimuthal distribution 0*. The inner error bars show the statistical errors, while the 
outer error bars represent the sum of the statistical and systematic errors added in quadrature. 
The solid lines show the results of the two-dimensional fit to the data (see text). The dashed 
lines show the expectations for s-channel helicity conservation (SCHC). 
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Figure 5: The three spin density matrix elements a) r^, b) rjij and c) Re [r^] for p meson 
photoproduction as a function of \t\ (full points) together with ZEUS measurements [3] (open 
points). The inner error bars show the statistical errors, while the outer error bars represent 
the sum of the statistical and systematic errors added in quadrature. The full lines show the 
predictions of the BFKL model and the dotted lines show the predictions of the two-gluon 
model. The dashed lines show the expectation from s-channel helicity conservation (SCHC). 
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